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In' the studies that have thus far appeared on the 
aerodynamics of ro tat ing-v/ing aircraft, rotor systems are 
investigated of such a character that rotation of the 
"blade a "bout its span axis, except for torsional deflection, 
is excluded from consideration. In the present report, 
with the aid- of the usual computation methods, a rotor is 
investigated the pitch of whose blades is capable of feeing 
controlled in such a manner that it varies linearly with 
the flapping angle. To test the effect of this linkage on 
the aircraft performance, the theory is applied to an illus 
trative example. 



X. PRELIMINARY REMARKS 
1. Control Linkage 



The aerodynamics of r o t at ing-w ing aircraft has been 
treated in a considerable number of published reports. 
These arc all concerned with wing systems for which the 
rotor blades, hinge-connected to the axis of rotation, 
possess tvo degrees of freedom; a motion of rotation of each 
blade in the plane of rotation, permitted by the swiveling 
hinge in the rotor disk (plane containing the hinge and the 
perpendicular to the rotor axis), and a flapping motion at 
right angles to the rotor disk permitted by the flapping 
hinge. With this arrangement, no deflection of the blade 
occurs about its span axis* The blade pitch, that is, the 
angle between the zero lift line of the blade section and 
the plane of rotation (fig. 3) is thus determined by the 
design and remains constant during the rotation of the 
wing system. 



* n Acrodynamik des Drehfluglers mit Blat twinkelriicks touerung. 
Luf tf ahrtf orschung, vol. 16, no. 7, Jul; 20, 1939, pp. 355- 
361. : 
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In recent times, a different kinematical system has 
proven itself practically applicable* With this system, 
a rotation of the wing about its span axis, in addition 
to flapping and rotating, is possible. This variation of 
the blade pitch is, however, positively coupled with the 
flapping angle in such a manner that the upward motion of 
the blade reduces the blade pitch. The process can best 
be understood by considering the design of the rotor sys- 
tem in detail. Such coupling has been practically applied 
to the Breguet helicopter (G-erman Patent No. 567 , 584/ 62b 
(1933)) and the Hafner gyroplane (reference l) . Since the 
latter is described in detail in the literature, the ar- 
rangement of its rotor system has been used as the basis 
for our computation. 

The essential characteristic of the design is the 
so-called IT spider ir shown in the sketch of figure 1. It 
consists of the spider axis and the spider arms which are 
connected with the rotor blades through hinges G- and 
lever arms H. The rotor blades are attached to the ring 
H through a Cardan hinge K, making possible the change 
in blade angle with flapping angle. In the Hafner gyro- 
plane, the flapping hinge is through suitable means placed 
in the axis of rotation. In figure 1, the more general 
case of a distance e "between the hinge K and the axis 
is assumed. The spider axis may move up and down in ring 
R and be inclined in any direction. The possibilities 
thus provided for the change in the blade pitch are of 
great importance for the controllability of the aircraft 
but will not be considered further here. tfe are inter- 
ested only in the relation between the flapping and blade 
angle for steady flight conditions with constant setting 
of the spider axis with respect to the ring R. 'This re- 
lation may be derived with the aid of figure 1. Let the 
flapping angle be denoted by 3, the blade pitch by g 
(see figs. 2 and 3); for p = 0, let d = d 0 » For an up- 
ward motion of the blade point A at which the lever am 
H is attached to the blade spar, moves up by the amount 
p(a-e). Since the position of the hinge G remains un- 
changed, the blade pitch is decreased by the amount 
d — e 

P • Between $ and |3 , we thus have the linear 

relation: 




For briefness, we set 
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« m --*Jl± ■ ■ (1) 

R h 

and denote it as the linkage ratio. We thus have 

4 = £ 0 - * R § • ■ ; (2) 

With regard to the order of magnitude of v% , the 
following nay be said. If the flapping hinge K is situ- 
ated in the axis of rotation dg = a/h: i.e., equal to 
cotangent of the angle formed by the "blade axis with the 
line joining the hinge 3 and the center of the rotor 
system. For reasons 'of symmetry, this angle will he 
chosen preferably about 60°. The linkage ratio will then 
have a value of dg = 0.577. If e is greater than 0, 
a smaller value of v R will he obtained; for e * a, 
d« = 0 and we obtain the usual system of the Cierva auto- 
giro. In general, e should be as small as possible, so 
that in practical designs it will lie between the values 
0.3 and 0.6. 

2. Statement of the Problem and Symbols Used 

The object of the following computation is to extend 
the results already obtained on the aerodynamics of rotat- 
ing-wing aircraft in steady, forward flight to the case 
where we have the above-described linkage between the blade 
and flapping angles. 

The papers on gyroplanes are all based on the original 
investigations of Glauert and Lock. Following upon the 
work of these men, a considerable number of further inves- 
tigations have been conducted of which chief mention is to 
be made of a paper by V/heatley (reference 2) in which the 
entire theory of the computation of the air forces of ro- 
tating-wing aircraft is built up anew and the old computa- 
tions arc considerably extended. 'There is also to be men- 
tioned a paper by Sissingh (reference 3) which follows 
closely the work of tfheatley but introduces further refine- 
ments, so that it is possible to estimate the importance 
of certain of the assumptions of Wheatley by a more accu- 
rate computation. In Sissingh f S paper will be found a de- 
tailed bibliography, so we shall here dispense with it. 



The computation given below largely follows the work 
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of Wheatley. We may therefore dispense with a detailed 
explanation of the assumptions and formulas which hold for 
all rotating-wing aircraft independent of any particular 
kinematics of the wing system. The investigation is re- 
stricted to the simplest case of untwisted rectangular 
"blades: i.e., the blade pitch $ 0 and the blade chord are 
constant* The "blade twist due to torsional moments about 
the "blade-span axis is by suitable constructive means 
( smallest possible distance between "blade center of gravity, 
shear center of blade cross section,, and center of pressure 
of blade section) held low enough so that it may he ne- 
glected. 

The following notation is used: 

R, rotor radius (distance of blade tip from rotation 
axis ) ; 

z, number of rotor blades; 

t , blade chord ; 

zt i . , . , 

a = — , solidity; 

17 H 

V, forward velocity; 

uu, angular velocity of the rotor blades; 
v/, induced velocity; 

a, angle of attack of the rotor "disk. 

Then 

V cos cc is the component of the forward velocity in 

the rotor disk; 

V sin a - w, component of the forward velocity normal 

to the rotor disk (fig* 2). 

Prom these, we further define the nondimens ional coeffi- 
cients : 

V cos a 

\i = . 

w E 

. V s i n a - w 

A = 

uj E 
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A, 


rotor 


w, 


rotor 


s, 


rotor 



For the purposes of our investigation, wo may with 
sufficient accuracy set S equal to the resultant air 
force, so that A = S cos aj 

Mdj applied moment of the rotor blades about the 

axis of rotation (for the autogiro ll& = 0) ; 

S ... 
lc s at , thrust coefficient; 

1 2 - 4 

~ p OJ R TT 

w 

c s — , lift-drag ratio; 
A 

Ihe above-defined magnitudes refer to the rotating- 
wing system as a whole. ?or the flow and lift relations 
at each blade element, the following symbols are used. 
(See fig. 3 . ) 

r, distance of a blade element from the axis of rota- 
tion; 

r 

* " E 

m, blade mass per unit length along blade span; 
E 

J = / m r 2 d r, mass moment of inertia of a blade 
' 3 about the poir. t . r ■ 0; 

H - I m g r d r, moment of the blade weight about 

B J the point r = 0; 

^ 0 

a r , angle between zero lift line and resultant flow 
at blade element at position r; 

e $ blade section lift coefficient; 
c T , constant mean value of slope of curve c a - f(a r ) 
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c w , constant mean value of drag coefficient c w ' of 
"blade section; 

Up 0)1, velocity at blade element parallel to rotor 
axis; 

Ui U)S t velocity at blade element in the rotor disk; 
u t 

cp = — , angle between direction of resultant flow and 
u p rotor d i s k ; 

, azimuth angle of "blades measured in plane of 
rotor (fig. S) 5 

ot"£ 4 c ! 

V = - — — , "blade mass constant; 

J 

M s , moment of air forces of a "blade about the flap- 
ping hinge ; 

B, factor to take account of the tip losses. 

The ratio |i serves in the computation as the inde- 
pendent variable. The object of the investigation is to 
determine the flapping angle |3 and the magnitudes k g , 
a, and €* 



II. THE AIRCRAFT IN STEADY FORWARD FLIGHT 
1. Induced Velocity -and Angle of Attack 



In determining the induced velocity and the angle of 
attack, the rotor system is considered as a lifting vortex 
of width 2R. Under the usual assumption that the induced 
velocity w and hence also the ratio A are constant for 
the entire blade—swept area, we have, according to Wheatley 

w = , uj S (3) 

4 7u 3 + X 3 

end 

X k s 

tan a = — + ■ ■ (4) 

k c and A arc for the present unknown magnitudes and 
will be determined below as functions of |ju (It is to 
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"be noted that between the thrust coefficient k s and the 
magnitude Crp employed by Who at ley there is the relation 
k s * 2C T ) 

2. Flow Velocities at Blade Element 

For the velocity at the "blade element in the plane of 
the rotor, wc obtain 



(5) 



u t = x + (i s in 

and for the velocity parallel to the rotor axis 

u-n U>E = X u> E - r — - - M> UJ E £ cos ^ 

x d B 
* UJ d t 

where Cl 1 is the derivative of the flapping angle g 
d t 

with respect to tine. The angle g (acute angle between 
the blade span axis and the rotor disk) is expressed in 
terns of the azinuth angle y "by a Fourier series in which 
terns only up to the second harnonic are retained: 

P = a 0 - a x cos \|/ - hi sin V - a 2 cos 2 V - b 2 sin 2 i r (?) 

•There is then obtained from (6), with \j/ = U3 t 

I / 1 \ i 

u = X + i p, a x + ( - p, a Q + x "b i + _ p a 2 jcos V 
P P \ 2 7 



+ (^~xa 1 + ipD 2 ^siny -f^~jjia 1 + 2xh 2 Jcos 2 ^ 

+ i jj> *b x ~ 2 x a a sin 2 ^ 
2 

+ ht a, cos 3t + - ^ ^2 sin 3 ^ ( 8 ) 

2 2 



The radial velocity in the direction of the "blade span 
axis need not be taken into account. Its effect on the 
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magnitudes to 'oe investigated will as usual "be neglected. 

3. Flapping Angle 

The flapping motion j3 is obtained by determining 
the Fourier coefficient a 0 , ai, bj, an b 2 « In most of 
the investigations on the aerodynamics of helicopters and 
aut ogiros , the second harmonics a s and b 2 , in contrast 
to the procedure of Wheatley, are dropped, "being consid- 
ered negligible. This is justified "by experience as long 
as v/g are mainly concerned v/ith the aerodynamic behavior 
of the rot at ing-wing system In its total effect: i.e., in 
the magnitudes S, A, tf, u), V. For the designer, it is 
important, however, in order to estimate the strength of 
the rotor blades, to know as accurately as possible the 
air forces at each blade element which forces are con- 
siderably affected by the second harmonic of the flapping 
motion. The coefficients a 2 and b 2 .arc therefore in 
the following retained. 

The differential equation of the flapping motion is 
obtained from the moment equilibrium equation about the 
flapping hinge. In this condition there enter the air 
forces, the inertia forces of the flapping blade, the cen- 
trifugal forces, and the weight • With the notation given 
above, there is obtained the known relation 



d s j 



+ 3 



(9) 



It is here as sumod that the distance e "between the flap- 
ping hinge and the axis of rotation is sere or small enough 
so that its effect may "be neglected. 

The moment of the air forces about the flapping hinge, 
under the assumption that the resultant velocity at the 
blade element is v/ith sufficient accuracy given "by u^.U)R, is 

35 

M s = J I p t u t 2 w 2 H 2 c a r d r 

J 2 

0 B 

M g = i p t uj 2 R 4 / u t 2 c a x d x (10) 

*6 
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To take account of the thrust losses at the "blade tip, we 
integrate, following Wheat ley, not up to 1 hut only 
up to BE, v/here 

t 

B = 1 - 



ETow 



and from figure 3 



2 R 



a r = i} + cp 



With the aid of (2), there is then obtained 



This holds, however, only for the advancing blade section. 
For the retreating "blade (it < vjj < 2 tt and 0 < r < - 
|ju r sin\|/), co r is as usual replaced "by -a r and it is 
assumed that the slope c a 1 of the lift curve has approxi- 
mately the same value as for the advancing blade* We set 

Mg = Kg + Kg 

v/here Mg 1 is the air-force moment that is obtained by 

substituting the angle a r valid for the forward moving 
"blade, and Kg is the correction term that takes the 
retreating motion into account. From (10), there is ob- 
tained 
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Ml 1 



= 0 for 0 < y < * 
— m sin 



- - t , tur ]<• c a ' JV, - 0.0) Mr + u, u,] *<*x 



*i — -f ^ * a a. + - - * &l>i-\ B* a t 



(12 
a, b) 



(15 c) 



+ 7 (f M 83 Ul + B * b * + /<2 B * bi ) 
With the aid of (7), there is obtained from (9) 



for n < y < 2tz 

Ifg it now developed into a Fourier ser- 
ies of which again only terms up to the 
second harmonic are retained: 
M 8 ■» -i- q tv? & c a ' {A 9 + A x cos y> + B x sin v 

+ A t cos2 y> + £ g sin2v) . . (13) 
For each of the harmonics, if the small effect 
of the retreating motion is taken into account 
only in the coefficients A 0 , A x , and B x , the 
following expressions are valid: 

1 r m ° 

^■5 t— ^ 

o 

. i n _MV_. 

^ ~i — i^TT 

J J etw'B'c,,' 

M 

2n 

1 Af| (cos v, sin v) 
71 \ 1 



1 Af? (cos y, sin v) 

+ « \ — *» 



1 f 3f J (cos 2 v, sin 2 y) 
n — dv 



(14 a— c) 



Substituting (5), (8), and (7) in (12) 
and then in (14), there is obtained, 
after carrying out the integrations 

^-"l^esr-^l^ + T** 



- 2?« b x - | 5' b x + y ai * s a,) , 
= 1 /i B* a x + j 2?« S, - ~ ^ * 3 0 O 



(15a-d) 



^QtattR* c a ' 



|y |a 0 -h 3 a a cos 2 v + 3 6, sin 2 y — ^j^ jj (16) 



With 



y = 



Q t R* C q 



By comparing the coefficients of the trigo- 
nometric functions of equations (13) and 
(16), the following relations are obtained. 



2 . 2M ff 

- a 0 H j~y 

y yJ co 2 



A 0 , 



A x = 0, 
B x = 0, 



6 



(17 a— e) 



— a 3 = A t , 

y 

V 

Iquations (17) in connection with (15) 
are now to be solved for the Fourier 
coefficients of the flapping motion. 
It is sufficient here to compute a 0 , 
a i» t> i up to terms of the order of 
magnitude M- and a 5 , b a up to terms 
of the order of p. 2 . Under this assump- 
tion a 3 and b s may be determined 
from equations (17 d,e) if, in the ex- 
pressions for A a and Ba, the follow- 
ing approximation values are substi- 
tuted for the coefficients a 0 , a x , b x 



a*' 
a x ' 



2/i 



^(l+tf, 



4 <*<> n 



(18 a— c) 



3 ^ B 

These relations are obtained from (17), 
if in (17a) for the computation of a 0 
only terms without \x are retained and 
in (17 b t c) for computing a lt b x only 
terms with \± l . With the aid of (18). 
there is then obtained from (17 d,e) 

B* 



72 + 9y^B* 



a % = ■ 



lyfi 2 B 



» (19a,b) 



18 (8 + y^ B*) 

Knowing a 8 and b a , the remaining co- 
efficients of the flapping motion may be 
determined from (17 a-c) up to the order 
of magnitude of p. 4 . Carrying out the 
computation, there is obtained 



l 

1 ^ b 1 v [ m \±- + 2 tvj - 2 ,i 

+ <P t - * [3 A f? 2 + -?-A^ + 4^ 0 B (^ + /^ a ; 
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with 0 l = (^-f-4^)( l + ^ , ) + 2 ^^ 



+ y^« a )-j^6 l (H-2V)]. 



4/i 



4. Thrust Coefficient 

For the total thrust of all of the 
t blades, there is obtained for one 
revolution the mean value 

2* Bit 

0 0 

For the thrust coefficient, there is obtained 

|« B 



= kS S Ui%Cadxdw - 



where 



2n 



Substituting the angle of attack a r and sub- 
dividing the range of Integration for advancing 
and retreating motion, there is obtained. 

2i B 



*s = ^ { J J [(*. - 0) Uit + utujdzdy, 



0 0 
2n — fisinxp 



(20) 



(21) 



With the aid of (5), (8), and (7) 9 there 
is obtained from (20) after evaluating 
the Integrals up to terms of the order 
of magnitude M- 

5. Drag-Lift Coefficient 
To determine the drag and power of 
the rotor a knowledge of the drag-lift 
ratio c is required, the latter being 
computed from the energy losses of the 
rotating-wing system. These losses 
arise from the induced velocity and from 
the drags at each blade element in ro- 
tating. The first-named portion of the 
loeses is obtained from considerations 
on the rotor system in its total effect, 
the corresponding formulas not being 
effected by the linkage between the 
pitch and flapping angles. The linkage 



ratio 



R 



does not enter into the 



second portion of the losses because a 
mean value of Cg independent of the 
angle of attack is used in determining 



the drag of a blade element . Wheat ley 1 s 
value of the lift-drag ratio can there- 
fore be taken over without modification, 
his expression being 

6. Balance of Moments about the 

Vertical Hinge 

All of the magnitudes thus far de- 
termined contain the nondlmensional ra- 
tio X whose dependence on \i is still 
to be determined* A relation between A 
and M- ie obtained from the condition 
of moment equilibrium about the vertical 
hinge. This condition is 

BR 

o o 

'iir R 

— ^j J jQtufor K-c v rdrd V = 0 



or 



4- J J* w<* c a <p x d x d y> 



-J* J* uf c w x dx d tp = 0. 



The first of the double Integrals 
represents the contribution of the lift 
to the moment, the second, the contribu- 
tion of the drag. The latter integral 
is to be Integrated from 0 to 1, not only 
to B because the drag is not decreased 
by the thrust drop at the blade tips. 
Taking account of the retreating motion 
of the blades, there is obtained after 
substituting the expressions for c a and cp 

2.i B 

Q 0) > /f 5 a Ca o o 
2f — n sin M> 



+ J JVo-<>*0) u, u t + u v *)xdxdy, 



2^ — n sin M> 
-2 J /[(*<>— 0*0)"p"<+ u v *]xdxd 



V 



(23) 



0 0 

2-1 i 2^ —ftslaxp 

— °fr J J* u t *xdxdtp-\-2^ J J ufxdxdy 
° 0 0 ° 0 0 

Performing the integrations and rearrang- 
ing the terms, there is finally obtained 

- § M B> o, 6, + I „« P (o,» + 1 «,» - o, a,) 



+ B«(o,»+6,«) +± M B> (0,6,-6,0.) + { #./.'B'6, 
- »„ [ J, fl'a, a, + 1 ff B» (a, 6, - 1 a, 6.) 
-ifi. Ma, o, + 6, ^J-^jl+^-i^) 



+ 



= 0 



: £ co 1 i? 5 7i a c a ' 



(24) 
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The practical computation of the ratio A as a func 
t ion of \i is somewhat difficult since (24) contains the 
coefficients of the flapping motion, which coefficients, 
according to (19), again depend on A and \x* If A is 
known, however, no difficulty is encountered in determin- 
ing (3 , a , k s , and G • 



III. ILLUSTRATIVE EXAMPLE 



The formulas set up in the preceding section will nov 
"be applied to a numerical example. The object of the com- 
putation is not so much to clarify the process of practice 
computations as to obtain a numerical estimate of the sig- 
nificance of the pitch and flapping angle linkage on the 
aire raf t pe rf ormanc e . 

As an example, there will "be chosen an autogiro of 
the usual present-day design. The initial values used as 
a basis for the computation and which approximately corre- 
spond in their order of magnitude to the design factors 
used for the C 30 type autogiro are the following: 



G- = 900 kg M d ■ 0 

R = 6.00 hi c a * = 5.6 

z = 3 c w = 0.014 

t = 0.28 m # 0 = 6° 

J = 2 5.9 mkg s 2 d R = 0.45 

B = 0.9 3 

With these values, there is obtained from (19) in 
connection with (24) the variation of the Fourier coeffi- 
cients of the flapping motion with |i shown in figure 4 
In figure 5, k s is plotted from (21) and in figure 6, 
the lift-drag ratio € from (22). For A = G, 



cos a 
m o 

the rotational speed n = — u) is obtained from 

TT 
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k = 

1 

2 P OJ 2 H. TT 




2 Gr 



fc s p tt E 4 cos a 



The result is plotted in figure 7 for G = 900 kg. The 
depender.ee of the forward speed on p>, which is obtained 
from the relation 

„ _ M> g J 

cos a 



is shown on figure 8. The maximum velocity of the aircraft 
will lie near the value \i = 0.4 (V = 208 kxa/h) . For this 
flight condition, figure 9 shows the corresponding flapping 
m o t i on . 

The coefficient a c , denoted as the coning angle, has 
the value 4.73°. To a flapping angle of this magnitude, 
there corresponds the blade pitch angle 

= a o = ^o m a o 

= 6 - 0.45x4.73 



m 3.87° 

The rotor thus operates on the average with a pitch of the 
usual order of magnitude of about 4 J . 

The important result brought out in figures 5 to 9 
is that the autogiro investigated having the blade pitch 
control described does not show any particular, unexpected 
properties but behaves entirely like a normal autogiro of 
the Cierva type. (Compare the computations of Wheatlcy, 
U.AtCnAi Reports ITos. 487 and 591.) Particularly note- 
worthy is the fact that the amplitude of the flapping mo- 
tion, i.e., chiefly the coefficient a % , is of the usual 
magnitude, so that the type of linkage described does not, 
as might have been expected, lead to a decrease in the 
flapping motion. 

All this leads to the final conclusion that the kine- 
matics of the blade-pitch control of the autogiro, while 
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offering no disadvantages, does not offer any particular 
advantages as far as "behavior in steady forward flight is 
concerned* It is true that in flight under gust condi- 
tions the linkage gives the rotor less sensitivity to the 
gust. The use of the linkage will hardly prove worth-- 
while, however, unless other advantages are gained through 
its application. The advantages of the spider construc- 
tion described in section I are to he found chiefly in the 
control possibilities of the r otat ing-wing system, "because 
of the fact that the spider axis may "bo displaced up and 
dov/n or tilted to the side in the ring H. The effect of 
the linkage is to he considered a by-product conditioned 
by the design which, however, does not impair the aerody- 



namic behavior of the aircraft, 
however, that caution must he empl 



It should be pointed out 
oyed in generalizing 



the results obtained from the computation example 



age ratios that 
inves t igat ed« 



considerably exc c c 



IV. SUMMARY 



The m 
of the air 
in that in 
dependence 
angle. Th 
autogiro w 
that for s 
on the aer 



ethod developed by Wheat ley for the computation 



d the value of 



to 



link- 
a 0.45 



forces of rotat ing-wing 
stead of a constant 

is assumed between thai 
e application of the 
ith a linkage ratio of 
toady, forward flight no 
©dynamic behavior of the 



aircraft is extended 
blac.e-pitch angle a linear 
angle and the flapping 
formulas obtained to an 

the usual magnitude shows 
particular effects arise 
aire raft ♦ 



Translation by S. Reiss, 
national Advisory Committee 
for Ae ronaut i c s . 
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Figure 1.- Autogiro rotor with linkage between the 
blade pitch angle and the flapping angle. 
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Figure 2.- Velocities and forces of the rotor system. 




Figure 3.- Velocities at the blade element. 
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Figs. 4, 5,6,7, 8, y 
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Figure Rotational speed as a 
function of [i for G = 
900 kg . 
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Figure 4.- Fourier coefficients of 

the flapping motion as 
a function of the p, * 
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Figure 5.- Thrust coefficient as 
a function of jjl . 



Figure 8.- Forward velocity as a 
function of p, for & s 
t?00kg , 
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Figure 9.- Flapping angle as a func- 
Figure 6*- Rotor irag-lift coef- tion Of the azimuth angle 

ficient as a function f or g - 900 kg and p,= 0.4 * 
of ul . 



